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I.  INTRODUCTION 

The  field  of  adaptive  array  processing  has  been  in  existence  for  the  past  two 
decades.  Foundations  of  the  field  were  laid  by  Howells  and  Applebaum  with  the 
introduction  of  the  adaptive  sidelobe  canceller  in  1965-1966  [1,2]  and  by  Widrow 
et.  al.  with  the  development  of  the  broadband  adaptive  antenna  in  1967  [3].  Also 
in  1967  and  the  year  following,  Capon  et.  al.  and  Lacoss  published  articles  detail¬ 
ing  their  work  on  the  signal  processing  of  data  from  a  large  aperture  seismic 
array  [4,5].  Since  then,  a  variety  of  different  adaptive  array  processing  techniques 
have  been  described  in  the  literature  [5-15].  During  that  period,  a  type  of  pro¬ 
cessing  called  constrained  array  processing  was  developed  [4,6,7,10,14].  One  of 
the  more  notable  papers  on  this  subject  is  the  one  by  Frost  [7]. 

In  this  paper,  we  introduce  a  new  linearly  constrained  adaptive  array  pro¬ 
cessing  technique  which  performs  the  frequency-dependent  array  weighting  using 
filters  possessing  both  poles  and  zeros.  It  is  shown  that  the  need  for  such  weight¬ 
ing  arises  when  the  relative  bandwidth  of  the  received  signals  is  large. 

In  section  II,  we  begin  by  deriving  the  optimal  array  weighting  required  to 
eliminate  m-1  interference  sources  incident  on  an  m-element  array  with  a  con¬ 
strained  look  direction.  The  form  of  the  optimal  weighting  is  later  used  as 
motivation  for  the  development  of  an  adaptive  array  processor  based  on  pole-zero 
filters.  Also  in  the  section,  we  analyze  the  effect  of  isotropic  thermal  noise  on  the 
optimal  weighting.  Concluding  this  section  is  a  brief  discussion  of  the  relation¬ 
ship  between  a  constrained  processor  and  a  Howells-Applebaum  processor. 
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In  section  HI,  the  issue  of  approximating  the  optimal  array  weighting  is 
addressed.  The  all-zero  and  all-pole  filter  models  are  presented.  Section  IV  intro¬ 
duces  an  alternate  implementation  of  a  conventional  Frost  array  processor  and 
analyzes  its  performance  as  a  function  of  the  number  of  taps.  Simulations  of  the 
conventional  Frost  array  processor  are  presented  for  the  purpose  of  comparing 
with  the  pole-zero  Frost  array  processor  developed  in  section  V. 

In  section  V,  problems  associated  with  the  adaptation  of  pole-zero  filters  are 
discussed.  It  is  shown  how  a  simple  reformulation  of  the  adaptive  filter  error 
function  eliminates  many  of  these  problems.  Two  undesirable  attributes  of  the 
error  reformulation  are  discussed  along  with  ways  of  overcoming  them.  Simula¬ 
tions  are  presented  which  show  that  the  new  pole-zero  processor  can  offer  a  sub¬ 
stantial  improvement  in  the  array  performance  compared  to  that  obtained  using 
the  conventional  all-zero  method. 


-  4  - 

0.  OPTIMAL  CONSTRAINED  ARRAY  PROCESSING 

Overview 

A  constrained  array  processor  such  as  the  one  proposed  by  Frost,  operates  in 
an  environment  where  a  single  desired  signal  is  incident  upon  an  array  from  a 
known  direction  and  several  interference  signals  are  incident  from  other  unknown 
directions.  The  objective  of  Frost’s  processor  is  to  minimize  the  response  of  the 
array  in  the  direction  of  the  interference  signals  while  leaving  the  response  in  the 
look  direction  unaltered  or  fixed  at  a  specified  value.  A  block  diagram  of  a  Frost 
array  with  a  unity  gain  look-direction  constraint  is  shown  in  Figure  1. 
Throughout  this  paper  Z-transform  notation  is  used  to  indicate  a  fixed  linear 
digital  filter.  Conversion  from  the  Z-domain  to  a  function  of  frequency  is  accom¬ 
plished  through  the  substitution  z—e^7,  where  u  is  the  radian  frequency  and  T 
is  the  sampling  period.  In  Fig.  1,  the  look  direction  is  assumed  to  be  perpendicu¬ 
lar  to  the  array  axis.  A  non-perpendicular  look-direction  can  be  accommodated 
by  adding  alignment  filters,  usually  in  the  form  of  steering  delays,  to  each  sensor 
to  cause  the  look-direction  signal  to  appear  in  time  coincidence  at  the  output  of 
each  delay.  Such  a  system  is  often  referred  to  as  a  signal- aligned  array  [17]. 
After  alignment,  each  signal  is  passed  through  a  linear  filter  which  serves  as  the 
frequency-dependent  array  weight.  The  output  of  all  the  filters  are  summed  to 
form  the  overall  array  output.  Because  the  look-direction  signal  appears  identi¬ 
cally  at  the  input  to  each  array  filter,  it  will  experience  a  response  through  the 
array  determined  by  the  sum  of  all  the  array  filters.  This  sum  is  constrained  to 
unity  in  order  to  maintain  an  undistorted  response  in  the  look  direction.  Subject 
to  the  unity-sum  constraint,  the  filters  are  adjusted  to  minimize  the  output 
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Figure  1.  Block  diagram  of  a  Frost  array  with  unity  gain  look-direction  constraint. 
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power  of  the  array  and  in  this  manner  eliminate  interference  signals  without  dis¬ 
torting  the  look-direction  signal. 

Optimal  array  weighting  --  Deterministic  formulation 

We  now  derive  the  optimal  frequency-dependent  array  weighting  required  to 
cancel  m-1  broadband  directional  interference  sources  incident  upon  an  m-sensor 
array  with  a  constrained  look-direction  response.  The  following  analysis  assumes 
ideal  propagation  conditions  and  negligible  thermal  noise.*  Let  the  Z-transform 
of  the  look-direction  signal  be  S(z)  and  the  transform  of  the  ltk  interference  sig¬ 
nal  be  Ri(z).  Since  the  look-direction  signal  appears  identically  at  all  sensors  of 
the  array,  maintaining  zero  distortion  in  the  look-direction  requires  thao 

s(0  £  WHO  =  S[z)  .  (1) 

t«*l 

Cancelling  S(z)  from  both  sides  of  eq.  (1)  leaves 

£  WHO  =  1  ,  (2) 

k—l 

the  look-direction  constraint.  For  simplicity,  assume  a  uniformly-spaced  linear 
array.  With  such  an  array,  each  interference  signal  undergoes  a  propagation 
delay,  At,  between  adjacent  sensors.  This  delay  is  determined  by  the  wave  pro¬ 
pagation  velocity  c,  the  sampliug  rate  T ,  the  intersensor  spacing  d,  and  the 
difference  in  angle  between  the  interference  direction  and  the  look  direction  fy, 

*  In  this  paper  any  signal  that  is  uncorrelated  from  sensor  to  sensor  is  referred  to  as  ther¬ 
mal  noise. 
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A  =  — — —sin^  . 
c-T  1 

Complete  cancellation  of  all  interference  signals  requires  that 

E  Wk(2)  z-{k-1^'  =  0  (3) 

i-t 


for  /= 1,2,  •  •  •  ,m-l.  Equations  (1)  and  (3)  can  be  written  in  matrix  form  as 


IT 


A(a)W(i)  = 


l  0 


where, 


S(z) 

RiU)  o 


i  l 

1  z~Al 


A(z)  = 


0 


v 


and 


W(») 


Wx(z) 

W2(z) 


(4) 


1 

z~(m~  1)A, 


5) 


(6) 


By  induction,  it  can  be  shown  that  the  solution  to  this  set  of  linear  equations  is 
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1 

5}  all  single  advances 
all  advance  pairs 

jjm-l.2(A!+  •••  +  A„_i) 

The  above  frequency-dependent  weighting  applied  to  a  linear  uniformly  spaced 
array  would  cause  all  of  the  interference  signals  to  be  perfectly  nulled  regardless 
of  their  bandwidth.  Notice  that  the  optimal  weighting  contains  a  set  of  m-1 
pole-like  resonances  common  to  all  sensors  of  the  array.  This  fact  will  be  used 
later  in  the  paper  when  developing  a  pole-zero  array  processor.  Also  note  that 
since  the  intersensor  delays  can  be  fractional  and  the  frequency  resonances  caused 
by  the  poles  are  infinite,  an  exact  rational  realization  of  this  weighting  is  not  pos¬ 
sible.  Instead,  one  must  be  content  with  a  rational  approximation  made  over  a 
selected  bandwidth.  The  issue  of  approximating  the  optimal  weighting  will  be 
addressed  further  in  the  next  section. 

To  gain  insight  into  the  frequency  response  of  the  optimal  weighting,  con¬ 
sider  a  single  inteference  source  incident  on  a  two-sensor  array.  In  this  case, 
m=2  and  the  optimal  weight  transfer  functions  are: 

»',(*)  =  and  IV2(.-)  =  .  (8) 

A  Z  L  —  Z 

Typical  values  of  A  fall  between  +  2  and  -2.  Figure  2  shows  one  period  of  the 
frequency  response  of  the  optimal  weights  found  by  evaluating  eqs.  (8)  with 
z=e,uT .  Normalized  frequency  is  defined  a s  uT  and  ranges  between  0  and  ic/ A. 
From  these  plots  it  is  evident  that  the  magnitude  of  the  optimal  weight  transfer 


Figure  2.  Frequency  response  of  optimal  array  weights  for  a  two-element  Frost  array  with 
single  interference  source  incident  on  the  array.  The  jammer  arrival  angle  and 
sampling  rate  determine  the  intersensor  delay  A.  (a)  Sensor  1  (b)  Sensor  2. 


-  10  - 


functions  are  equal  and  dependent  upon  frequency.  Also  it  can  be  seen  that  the 
phase  responses  are  linear  and  of  opposite  slopes.  The  infinite  resonance  at  zero 
frequency  indicates  that  it  is  physically  impossible  to  maintain  both  unity  d.c. 
gain  in  the  look  direction  and  a  d.c.  null  in  the  interference  direction.  Notice  also 
that  reducing  the  value  of  the  intersensor  delay  A  causes  an  expansion  of  the  fre¬ 
quency  scale  of  this  plot. 


Effect  of  thermal  noise  on  the  deterministic  solution 

When  considering  the  effect  of  thermal  noise  on  the  array  output,  it  turns 
out  that  in  some  situations  the  weighting  given  by  eq.  (7)  will  perfectly  null  the 
directional  interference  sources  but  at  the  expense  of  exagerating  the  thermal 
noise  contribution  to  the  array  output.  Thus,  while  the  signal-to-interference 
ratio  is  dramatically  improved,  there  could  be  a  substantial  degradation  in  the 
signal-to-thermal-noise  ratio.  Such  an  effect  is  easily  demonstrated  by  calculating 
the  thermal  noise  power  at  the  array  output  as  a  function  of  the  intersensor  delay 
A.  The  power  spectrum  of  the  array  output  due  to  thermal  noise  is  the  norm 
square  of  the  weight  vector  times  the  noise  power, 

1  y.l-H2  =  -y-iw(OI2  •  W 


For  the  two-sensor  case  presented  above,  eq.  (9)  becomes 


|V.(*)I*  = 


JV. 


(io) 


Eq.  (10)  is  plotted  in  Fig.  3  as  a  function  of  the  intersensor  delay  A  and  with 


N. 


-=.001.  Notice  that  as  A  becomes  small,  the  thermal  noise  contribution  to 


INTERSENSOR  DELAY 

Figure  3.  Broadband  thermal  noise  power  at  array  output  as  a  function  of  intersensor  dela; 

Array  weights  chosen  to  completely  null  directional  interference  without  regard 
to  thermal  noise. 


the  array  output  becomes  large.  As  A  becomes  large,  the  thermal  noise  power 
appearing  at  the  array  output  drops  to  a  level  equal  to  the  amount  of  noise 
added  to  each  sensor  (the  thermal  noise  floor).  Observation  of  the  converged 
antenna  pattern  explains  this  phenomena.  In  order  to  place  a  null  near  the  look 
direction,  the  array  antenna  pattern  billows  out  in  other  directions.  The 
increased  gain  in  non-look  directions  causes  the  thermal  noise  at  the  array  output 
to  become  large.  Thus,  when  the  interference  direction  is  close  to  the  look  direc¬ 
tion,  the  effect  of  a  small  amount  of  thermal  noise  becomes  important  and  must 
be  taken  into  account. 

Optimal  array  weighting  --  Wiener  formulation 

To  overcome  the  thermal  noise  problem,  a  compromise  must  be  made 
between  interference  nulling  and  thermal  noise  suppression.  The  Wiener  solution 
accomplishes  such  a  compromise  by  minimizing  the  overall  array  output  power 
subject  to  a  constraint.  Using  Lagrange  multipliers,  it  can  be  shown  that  the 
frequency-dependent  array  weights  which  minimize  the  mean  square  of  the  array 
output  subject  to  a  unity-gain  look-direction  constraint  satisfy  the  following 
equations  [7j: 

1)  (l 

1  1 

r„(*)-wu)  =  ;  /(i  i  i)r„'  ;  ,  (ii) 

l)  u 


where  R„(’)  is  the  m/m  spectral  covariance  matrix  of  the  array  signals  after 
alignment  and  W(;)  is  the  mxl  vector  of  optimal  frequency-dependent  weights. 
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To  connect  this  set  of  equations  to  the  deterministic  set  given  in  (4),  multiply 
both  sides  of  (4)  by  Ar(r') 

S(z)' 

0 

Ar(*-»)A(*)-W(*)  =  Ar(2-»)-  (12) 

.  0  , 

Interpreting  |S(z)|2  and  |/?/(z)|2  as  the  power  spectra  of  stationary  stochastic 
signals,  the  matrix  Ar(z_1)-A(z)  can  be  viewed  as  the  spectral  covariance  matrix 
associated  with  the  array  input  vector 

AV‘)-A(r)  =  R  „{z)  .  (13) 

Multiplying  out  the  right  hand  side  of  eq.  (12)  using  the  definition  of  A(z)  given 
in  (5),  it  can  be  shown  that 


It  is  also  easy  to  see  from  the  definition  of  A(z)  that 

(  10  .  .  .  Oj  A(z)  =  5(z)  [  1  1  .  .  .  1  ]  .  (15) 

Post-multiplying  both  sides  of  eq.  (15)  by  >t"1(z)  then  taking  the  norm  square  of 
the  result  gives 
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(l 

1 


[  1  1  ...  1  ]-A-*(z)- A-T(r') 


UJ 


|S(*)|: 


Combining  eqs.  (12,13,14,16)  gives 


1 

i 

1 

l 

R„(z)-W(z)  = 

• 

/(i  i  i) 

• 

.1 

i . 

(16) 


(ID 


the  Wiener  equations.  Thus  with  no  thermal  noise  present,  the  deterministic  for¬ 
mulation  given  in  eq.  (4)  is  equivalent  to  the  Wiener  formulation  given  in  eq. 
(11)- 


Wiener  solution  in  the  presence  of  thermal  noise 

Finding  a  closed  form  solution  to  eq.  (11)  in  the  presence  of  thermal  noise  is 
a  difficult,  if  not  impossible  task.  To  gain  insight  into  the  array  operation  in  the 

N, 

presence  of  noise,  an  exact  solution  can  be  found  for  a  two-sensor  array.  Let - 

2 

be  the  power  of  the  thermal  noise  relative  to  the  power  of  the  directional  interfer¬ 
ence.  This  relative  power  level  could  be  a  function  of  z  (frequency),  but  for  the 
following  analysis  it  is  assumed  to  be  constant.  In  this  case,  the  solution  to  eq. 
(11)  is 
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-(1+  ~)2'A  +  Z' 2A 

Wx(z)  =  - - - - - — 

1  -  (2+  N,)z~ A  +  z’2A 

and 


(17a) 


W^2(z) 


A.  A 

1  -  (1+  -  f)z -A 
1  -  (2+  Ne)z-A  +  2-2A 


(17b) 


Note  that  with  thermal  noise  present,  the  optimal  array  weighting  is  still 
represented  by  the  ratio  of  two  polynomials  as  it  was  in  the  deterministic  case. 
Also  note  that  when  No—0,  pole-zero  cancellation  occurs  and  eqs.  (17)  reduce  to 
eqs.  (8). 


For  the  two-sensor  case,  the  contribution  of  the  thermal  noise  to  the  array 
output  can  be  computed  using  eqs.  (9)  and  (17)  and  is 


ir.wi1  =  n. 


e+T'-"4 


-zA  +  (2+  N ,)  -  z-A 


(18) 


The  array  output  due  to  the  directional  interference  is 


IWI!  = 


1^,(2)^+  M^2(2)| 2 
1  +  2-A 


=  (  — )2 
v  2  7 


~2 A  +  (2+  Nt)  -  2~A 


(19) 


The  array  output  due  to  both  the  direction  interference  and  the  thermal  noise  is 
found  by  adding  eqs.  (18)  and  (19)  to  give 

iv.w(2)is  =  m<*)iJ+  mwi* 


j 
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*.  +  A)’ 

“  _,-A+  (2+NJ~*A  '  (20) 

To  determine  the  performance  of  the  optimal  two-sensor  array  for  a  wideband 
signal,  eqs.  ( 18-20)  can  be  integrated  over  the  bandwidth  of  the  signal. 

Formulas  (18)  and  (19)  were  numerically  integrated  over  the  frequency  range 
ff/4  to  3jt/4  and  are  plotted  in  Fig.  4  as  a  function  of  the  intersensor  delay. 
Compare  this  graph  to  the  one  presented  in  Fig.  3.  Notice  that  for  large  A,  the 
thermal  noise  dominates  the  array  output.  As  A  becomes  small,  interference  nul¬ 
ling  is  traded  off  for  thermal  noise  suppression  and  the  interference  signal  power 
dominates  the  array  output.  Thus  when  using  the  Wiener  weighting,  directional 
interference  nulling  is  compromised  for  the  sake  of  thermal  noise  suppression. 

Relationship  between  a  Frost  array  and  a  Howells-Applebaum  array 

The  Howells-Applebaum  array  processor  is  different  from  a  linearly  con¬ 
strained  processor  in  that  it  attempts  to  maximize  the  SNR  of  the  received  sig¬ 
nal.  In  doing  so,  it  allows  the  desired  signal  to  undergo  some  distortion  as  it 
passes  through  the  array.  In  this  sense,  the  Howells-Applebaum  processor  can  be 
thought  of  as  having  a  soft  constraint  in  the  look  direction  as  opposed  to  the 
hard  constraint  used  in  a  Frost  processor. 

The  converged  weight  vector  of  a  Howells-Applebaum  array  satisfies  the  fol¬ 
lowing  equation  [11], 

|  I  +  1-W(,)  =  .  (85) 

where  v  is  the  quiescent  steering  vector  and  7  is  a  constant  determined  by  the 
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dynamic  range  of  the  processor.  In  this  paper  we  have  assumed  the  array  to  be 
signal  aligned,  so, 


The  solution  to  eqs.  (25)  turns  out  to  be  similar  to  the  solution  expressed  by  eqs. 
(11);  in  the  limit  as  7  goes  to  infinity,  the  Howells- Applebaum  weight  vector  and 
the  Frost  weight  vector  differ  only  by  a  scalar  transfer  function.  This  similarity 
indicates  that  a  Howells-Applebaum  type  of  array  processor  could  also  benefit 
from  the  use  of  pole-zero  filters.  A  soft  constrained  array  processor  using  pole- 
zero  filters  is  not  investigated  in  this  paper  but  is  left  as  a  topic  for  future 
research. 
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ffl.  APPROXIMATION  OF  THE  OPTIMAL  ARRAY  WEIGHTING 

Role  of  the  relative  bandwidth 

As  was  noted  in  the  previous  section,  the  optimal  frequency-dependent  array 
weighting  is  not  perfectly  realizable.  An  approximation  of  the  optimal  weighting 
must  be  made  over  the  array  operating  bandwidth.  An  important  parameter  per¬ 
taining  to  the  difficulty  of  the  approximation  is  the  relative  bandwidth  of  the 
array  signals.  Relative  bandwidth  is  a  dimensionless  quantity  defined  as  the  ratio 
between  the  absolute  bandwidth  of  the  received  signals  and  their  center  fre¬ 
quency.  For  acoustic,  sonar,  and  seismic  arrays,  the  received  signals  are  at 
baseband  and  the  relative  bandwidth  is  large.  For  electromagnetic  arrays,  the 
received  signals  are  at  radio  frequencies  and  the  relative  bandwidth  is  much 
smaller.  In  any  case,  the  goal  of  the  adaptive  array  is  to  form  an  approximation 
of  the  optimal  weighting  derived  in  section  II  over  the  operating  bandwidth  of 
the  array.  Simulations  depicting  the  array  performance  as  a  function  of 
bandwidth  were  presented  by  Rogers  and  Compton  in  (18).  Also,  an  analysis  of 
the  nulling  performance  as  a  function  of  bandwidth  is  given  by  Mayv  an  et.al.  in 
[19].  Our  paper  is  concerned  not  so  much  with  the  array  performance  as  a  func¬ 
tion  of  bandwidth,  but  with  the  improvement  of  the  array  performance  at  a  fixed 
wide  bandwidth. 

Approximation  model 

Fig.  5a  shows  the  tap-delay- line  filter  model  commonly  used  to  perform  the 
array  weighting  in  conventional  adaptive  array  processors.  It  consists  of  a  series 
of  delays  whose  outputs  are  linearly  combined  through  the  multipliers  &,•  to  form 


the  filter  output.  The  transfer  function  of  this  filter  is, 


H(z)  =  b0  +  6,z  1  +  •  ■  ■  +  bnz  n  .  (22a) 


The  roots  of  the  H(z)  polynomial  are  called  the  zeros  of  the  filter. 

Fig.  5b  shows  the  direct  form  realization  of  an  all-pole  filter.  The  filter  out¬ 
put  is  delayed  and  fed  back  through  the  multipliers  a,-  and  summed  with  the 
filter  input.  The  transfer  function  of  this  filter  is, 


H(z) 


_ 1 _ 

1  +  djz'1  +  •  •  •  -f  a„z~B 


(22b) 


The  roots  of  the  denominator  of  H(z)  are  called  the  filter  poles.  A  cascade  of  the 
all-zero  and  the  all-pole  filter  generates  one  type  of  pole-zero  filter  structure.  In 
this  paper,  we  consider  the  use  of  this  type  of  pole-zero  filter  as  the  model  for 
approximating  the  optimal  array  weighting. 


Independence  of  array  weighting  on  frequency  shift 

It  will  now  be  shown  that  the  difficulty  in  approximating  the  optimal  array 
weighting  over  a  specified  bandwidth  is  independent  of  frequency  shift.  This 
implies  that  down-conversion  to  an  intermediate  frequency  (IF)  or  to  baseband, 
will  not  affect  the  performance  of  the  array  processor.*  This  is  most  easily  proven 
by  assuming  the  array  signals  to  be  analytic  and  the  array  filter  coefficients  to  be 
complex  (as  opposed  to  real).  An  analytic  signal  is  a  complex  valued  signal 
whose  negative  frequency  image  has  been  removed  by  combining  the  original  sig¬ 
nal  with  its  Hilbert  transform.  Recovery  of  the  original  signal  is  accomplished  by 
*This  does  not  consider  hardware  implementation  issues. 


Fig.  5b.  Direct  form  realization  of  an  all-pole  digital  filter. 
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taking  the  real  part  of  the  analytic  signal.  For  more  information  concerning  the 
analytic  representation  of  array  signals,  see  appendix  B  of  reference  [17]. 

Let  the  approximation  of  an  optimal  array  weight  be  given  by 


E 

m*)  =  - -  .  (23) 

E 

k—0 


where  «0=1.  Note  that  since  the  ak s  and  bk' s  are  complex  valued,  the  poles 
and  zeros  of  this  filter  do  not  necessarily  occur  in  complex  conjugate  pairs.  As  a 
function  of  frequency,  eq.  (23)  can  be  written  as 


W{u) 


(24) 


where  the  shifted  frequency  range  is  u;1<w<w2.  A  frequency  shifted  version  of 
this  transfer  function  is 


W(oj  +  uc) 


£  bke~i{lJ)+^)Tk 
£  ake-]'{u+u‘)Tk 


(25) 


where  the  pertinent  frequency  range  is  Wj -wc  <w<u2-wc.  Eq.  (25)  can  be  easily 
rewritten  as 


W{ijj  +  ue) 


£  bkeiw<Tke->uTi 
£  ake]U‘  Tk  e~ 3 uTk 


(26) 


Letting  VF' (z)  be  the  frequency  shifted  response,  eq.  (23)  becomes 
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W‘{2) 


S  V  *'k 
E  «*' 


(27) 


where 

6*'  =  bke~ju'n  and  ak>  =  «*«-*"• n  .  (28) 

Using  the  transformation  given  in  eq.  (28),  we  see  that  VU(r)  and  W’(z)  are 
equivalent  when  evaluated  over  their  respective  frequency  ranges.  This  proves 
that  approximation  of  the  optimal  weighting  over  a  specified  bandwidth  is 
independent  of  frequency  shift  and  implies  that  an  investigation  of  the  bandwidth 
characteristics  of  an  array  need  not  be  concerned  with  the  intermediate  frequency 


1 


of  the  receiver. 
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IV.  CONVENTIONAL  FROST  ARRAY  PROCESSING 

Alternate  implementation 


To  simplify  the  adaptivity  issue,  we  introduce  the  alternate  constrained 
array  processor  shown  in  Fig.  6.  This  structure  was  first  discussed  by  Applebaum 
and  Chapman  [13]  and  more  recently  shown  to  be  equivalent  to  a  Frost  array  by 
Griffiths  and  Jim  [16].  The  basic  idea  of  this  realization  is  to  use  a  preprocessor 
to  remove  the  look-direction  signal  from  the  adaptive  process  x-input  making  it 
impossible  for  the  adaptive  process  to  cancel  signal  components  from  the  array 
output  by  subtraction.  Interference  signals  appear  in  both  the  x  and  d  inputs 
and  can  therefore  be  cancelled  from  the  output.  Such  a  structure  maintains  the 
unity  gain  look-direction  constraint  without  using  a  constrained  adaptive  algo¬ 
rithm.  The  pole-zero  array  processor  developed  in  the  next  section  is  based  on 
this  alternate  implementation  of  Frost’s  processor.  It  should  be  noted  that 
several  different  types  of  preprocessors  were  discussed  in  [16].  The  type  shown  in 
Fig.  6  is  conceptually  the  simplest  but  is  not  necessarily  the  best  of  these. 

Transformation  between  the  two  array  processors  is  readily  seen  to  be 


^l(-) 

1 

1 

r  *,(*) ' 

\V2  (z) 

1 

-1  1  0 

Bfr) 

■ 

_  J_ 

• 

-1  . 

m 

• 

•  • 

0  .  1 

W^z\ 

1 

-1 . 

(29) 


For  the  m=2  case,  this  set  of  equations  reduces  to 


Wi(z)  =  1/2-  B(z) 
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W2(z)  =  1/2+  B(z)  .  (30) 

Solving  eqs.  (3)  for  B(z)  gives 

B(z)  =  1/2 1 W-  Wt[z) )  .  (31) 

The  optimal  transfer  function  given  in  eq.  (8)  thus  becomes 

B(z)  =  1/^-  ~  -  (32) 

1  -  z~a 

From  eqs.  (29)  and  (32)  notice  that  optimal  weighting  for  the  alternate  Frost 
array  processor  still  contains  a  set  of  poles  common  to  all  array  weights. 

Simulations 

We  now  present  simulations  of  a  two-element  Frost  array  processor.  In  sec¬ 
tion  V  of  this  paper,  these  simulations  will  be  compared  to  a  similar  set  of  simu¬ 
lations  of  a  pole-zero  Frost  array  processor. 

The  conditions  of  the  simulation  are  as  follows:  The  interference  signal  had  a 
normalized  center  frequency  of  7r/2,  a  100 %  relative  bandwidth,  and  a  bearing 
angle  45  degrees  off  the  look  direction.  There  was  no  thermal  noise.  The  array 
consisted  of  two  elements  spaced  one  quarter  wavelength  apart  at  the  signal’s 
center  frequency.  Such  a  configuration  corresponds  to  an  intersensor  delay  equal 
to  \/2/2.  The  array  filters  each  had  5  taps  (4  zeros)  and  the  tap  spacing  was  1. 

Figure  7  compares  the  frequency  response  of  the  converged  Frost  array 
filters  with  the  frequency  response  of  the  optimal  array  weighting.  The  cross- 
hatched  dark  lines  indicate  the  optimal  frequency  response  plotted  over  the 
bandwidth  of  the  interference.  The  continuous  curve  represents  the  frequency 
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Frequency  response  of  optimal  array  weights  (cross-hatched  marks) 
compared  to  5  weight  approximation  from  a  Frost  array. 

(a)  Sensor  1  (b)  Sensor  2. 
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response  of  the  4-zero  Frost  array  filter.  It  is  interesting  to  note  that  if  the  rela¬ 
tive  bandwidth  were  small  enough,  the  optimal  frequency  response  would  need  be 
matched  in  magnitude  and  phase  at  only  a  single  frequency.  Such  a  response 
could  be  realized  using  one  complex  weight  as  is  the  case  in  most  narrowband 
adaptive  array  systems.  In  this  paper,  however,  we  are  concerned  with  perfor¬ 
mance  over  a  wide  bandwidth  which  requires  the  use  of  a  digital  filter,  not  a  sin¬ 
gle  complex  weight  to  perform  the  array  weighting. 

Figure  8(a)  shows  the  broadband  antenna  pattern  formed  by  averaging  the 
array  sensitivity  over  the  frequency  range  of  the  interference  signal.  Figure  8(b) 
shows  the  frequency  response  of  the  array  in  the  look  direction.  Because  the 
look-direction  constraint  was  set  to  unity,  this  response  is  flat  over  all  frequen¬ 
cies.  Figure  8(c)  shows  the  response  in  the  direction  of  the  interference  signal. 
At  certain  frequencies,  the  attenuation  is  as  much  as  50  dB,  but  when  averaged 
over  the  entire  interference  signal  bandwidth,  the  attenuation  is  only  20  dB. 
This  simulation  thus  shows  that  the  Frost  array  is  working  fairly  well.  In  section 
V,  we  will  show  that  this  performance  can  be  dramatically  improved  using  a 
pole-zero  version  of  F rost’s  processor. 

Figure  9  shows  the  nulling  capability  of  the  Frost  array  processor  as  a  func¬ 
tion  of  the  intersensor  delay  or  equivalently  the  interference  bearing  angle.  The 
purpose  of  this  experiment  is  to  show  the  array  performance  in  the  presence  of 
thermal  noise  and  for  different  interference  bearing  angles.  The  condition  of  this 
simulation  were  the  same  as  those  in  the  previous  experiment,  only  -30dB  of 
thermal  noise  was  added  The  interference  bearing  angle  was  scanned  from  -90° 
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Figure  9.  Broadband  thermal  noise  power  and  interference  power  at  output  of  a 
two-element  Frost  array  plotted  as  a  function  of  intersensor  delay. 
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to  +90°  corresponding  to  an  intersensor  delay  ranging  between  -1  and  +  1.  At 
each  angle,  the  array  was  allowed  to  converge  and  that  the  individual  contribution 
to  the  array  output  due  to  thermal  noise  and  to  the  directional  interference  was 
evaluated.  Notice  that  the  directional  interference  dominates  the  array  output 
for  all  angles  unlike  the  performance  of  the  optimal  array  shown  in  Fig.  4.  This 
occurs  because  there  are  not  enough  weights  in  the  Frost  array  filters  to  perform 
the  proper  array  weighting  over  the  entire  interference  signal  bandwidth.  Again 
this  simulation  will  be  compared  in  section  V  to  a  similar  simulation  of  a  similar 
simulation  of  a  pole-zero  Frost  array  processor. 

Performance  as  a  function  of  the  number  of  taps 

Intuitively,  it  is  apparent  that  as  more  weights  are  added  to  the  array  filters, 
approximation  of  the  optimal  weight  frequency  response  will  become  more  accu¬ 
rate  and  the  null  depths  will  become  deeper.  It  is  instructive  to  determine  the 
performance  of  a  Frost  array  as  a  function  of  the  number  of  taps.  An  expression 
for  the  interference  null  depth  for  an  arbitrary  number  of  taps  can  be  derived 
when  the  interference  source  is  white  and  incident  on  the  array  from  an  angle 
corresponding  to  an  intersensor  delay  of  one.  In  the  configuration  of  the  previous 
experiment,  this  occurs  when  the  bearing  angle  is  90  degrees  off  broadside. 

An  exact  expression  for  the  null  depth  is  found  by  solving  the  Wiener  equa¬ 
tions  for  an  arbitrary  dimension.  The  Wiener  solution  of  the  array  shown  in  Fig. 

6  is  equivalent  to  the  Wiener  solution  of  a  Frost  array  [16]  but  is  somewhat  easier 
to  evaluate,  Determination  of  the  Wiener  equations  for  the  equivalent  Frost 
array  proceeds  as  follows:  Since  the  intersensor  delay  is  one,  the  Z-transform  of 
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the  signal  at  point  x  in  Fig.  6  is 

X(z)  =  R(z)-z~lR(z) 

=  ,  (33) 

where  R(z)  is  the  transform  of  the  interference.  The  power  spectrum  of  this  sig¬ 
nal  is 


|X(*)I’  =  |fl(.-)|!|-=+l  +  2-z'l  .  (34) 

Interpreting  \X{z)\“1  and  | /?(■?)) 2  as  power  spectra  of  stationary  stochastic  sig¬ 
nals  and  assuming  a  white  interference  source,  eq.  (34)  becomes, 

|X(z)|2  =  -z+1  +  2  -  ar1  .  (35) 

The  nzn  covariance  matrix  of  a  signal  with  a  power  spectrum  given  in  eq.  (35)  is 
a  tridiagonal  matrix  of  the  form 


2  -1 

-12.  0 
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0  -i 

-1  2  J 


(36) 


Similarly,  the  cross-correlation  vector  is 


P  = 


(37) 


The  Wiener  equations  are 
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R„  B  =  P 


(38) 


The  special  tridiagonal  structure  of  these  equations  allows  them  to  be  solved  for 
an  arbitrary  dimension  using  mathematical  induction.  The  result  is 


B  = 


n+  1 


l/2(n-l)) 
n-1 
n-2 


1 


(30) 


Using  eq.  (29),  this  array  weighting  can  be  converted  into  the  equivalent  Frost 
array  weighting; 


W,= 


n+  1 


r  1 

n 

-(n-1) 

n-1 

-(n-2) 

and  W2  =  — — 

n-2 

* 

1  n+  1 

.  -1  , 

,  1  . 

(40) 


The  array  response  in  the  direction  of  the  interference  source  is 


^.(*)+  *-'W2(:\ 


(41) 


Plugging  eq s.  (40)  into  eq.  (41)  gives 

W  =  ~~(l  +  +  •••+*-+I)  •  (42) 

n+  l 

The  power  gain  of  the  array  to  a  white  interference  source  is  found  by  taking  the 
sum  of  the  square  of  the  coefficients  of  polynomial  (42)  resulting  in 
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(  1+  1+  1+ 


+  1) 


1 

n+1  ' 

Eq.  (43)  indicates  that  the  array  power  gain  for  a 
inversely  proportional  to  the  number  of  taps  in  the 
Thus,  a  large  number  of  weights  is  required  in  order 
deep  wideband  null. 


(43) 

white  interference  source  is 
two-element  adaptive  array, 
for  the  Frost  array  to  form  a 
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V.  POLE-ZERO  FROST  ARRAY  PROCESSING 

Motivation 

Aside  from  adding  more  taps  to  the  array  Biters,  another  way  to  improve  the 
bandwidth  characteristics  of  the  array  is  to  use  array  filters  possessing  both  poles 
and  zeros.  In  Fig.  2,  we  see  that  the  frequency  response  of  the  optimal  array 
weighting  contains  infinite  resonances  at  the  frequencies  0  and  ?r/A.  Such  reso¬ 
nances  are  better  approximated  with  filters  having  both  poles  and  zeros  than 
with  zeros  alone.  In  this  section,  we  introduce  an  adaptive  array  processing  algo¬ 
rithm  which  adjusts  array  filters  possessing  both  poles  and  zeros. 

A  pole-zero  adaptive  filter 

Development  of  an  adaptive  algorithm'  to  adjust  pole-zero  filters  has  been 
the  subject  of  much  research.  One  of  the  more  popular  adaptive  pole-zero  filter¬ 
ing  algorithms  to  appear  in  the  literature  is  the  one  by  White  [28].  This  algo¬ 
rithm  uses  steepest  descent  to  find  a  minimum  of  the  mean  square  of  the  filter's 
error  function.  Another  popular  algorithm  is  commonly  called  recursive  max¬ 
imum  liklihood  (RML)  [23]  and  uses  a  Newton-type  method  to  find  a  minimum 
of  the  error  function.  Neither  of  these  algorithms  has  been  entirely  successful 
when  applied  to  problems  involving  real  data.  The  reasons  for  this  limited  suc¬ 
cess  are  manyfold.  First  and  foremost  is  that  the  mean  square  error  (MSE)  is  not 
a  quadratic  function  of  the  filter’s  coefficients.  This  results  in  an  MSE  function 
which  has  multiple  local  optima  in  addition  to  a  global  optimum.  Furthermore, 
the  MSE  becomes  infinite  when  the  filter  coefficients  are  adjusted  such  that  the 
poles  of  the  filter  lie  outside  the  unit  circle. 
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Some  insight  into  the  above  problems  is  provided  by  examining  the  simple 
single-pole  adaptive  filter  shown  in  Fig.  10.  Let  the  inputs  Xj  and  d;  be  station¬ 
ary  stochastic  signals.  The  MSE  of  this  system  can  be  found  in  the  following 
manner:  The  filter  output  is, 

Vj  =  xj  +  “iVi-i 

=  xj  +  (Xj-i  +  dil tj-2) 

=  xj  +  a,[ary_,+  al(x,_2  +  a,yJ_3)] 

=  etc.  (44) 


Accordingly, 


For  stability, 


The  error  c;  is 


Vj  =  Xj  +  a,iy_,  +  a?Xj  _2  +  •  •  • 

00  . 

=  E  a\Xj.i  ■ 

1=0 


l«i!  <  1 


00  / 

£y  =  di-Vj  =  di-  E  «i *j-i  ■ 

l  **0 


Thus  the  MSE  is  given  by 


(45) 


(46) 


(47) 


MSE  =  EVJ]  =  E{d/}  -  g  +  g 

t« 0  1=0 


g  £Iai+"  H8 


Defining  the  following  auto-correlation  and  cross-correlation  functions, 
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<*„(0  =  E\xix}+i\  I-*9) 

=  E\ Xj-di+t)  ,  (50) 

allows  eq.  (18)  to  be  written  as 

OO  00  oo 

MSE  =  E\df\-  2  £<.'*,<(/)  +  S  £  a'+ra^(/-m)  .  (51) 

/-0  1—0  m-0 

It  is  clear  from  eq.  (51)  that  the  MSE  is  not  a  quadratic  function  of  the  feedback 
weight  a,,  but  is  instead  a  very  complicated  function  that  could  have  local 
optima  and  be  very  difficult  to  search.  In  fact,  in  an  adaptive  filter  with  poles,  a 
general  analytical  solution  for  the  global  minimum  does  not  exist.  Also,  equation 
(51)  is  valid  only  as  long  as  inequality  (46)  is  satisfied;  the  MSE  goes  to  infinity 
when  (46)  is  not  satisfied.  Such  problems  do  not  arise  in  an  all-zero  adaptive 
filter  because  the  mean  square  error  function  is  a  quadratic  function  of  the 
weights  allowing  the  global  solution  to  be  found  quite  simply  by  solving  a  set  of 
linear  equations. 

Equation- error  approach 

An  adaptive  pole-zero  array  processor  could  use  the  RML  algorithm  to 
adjust  pole-zero  filters  at  each  sensor  of  the  array  depicted  in  Fig.  6  in  order  to 
minimize  the  array  output.  This  approach  might  work  for  some  applications  but 
turns  out  not  to  be  very  robust. 

Our  approach  to  adaptive  pole-zero  filtering  bypasses  many  of  the  above 
problems  usually  associated  with  pole-zero  adaptive  filters  by  reformulating  the 
filter’s  error  function.  This  is  done  in  a  manner  which  allows  the  polynomial 
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associated  with  the  poles  of  the  adaptive  filter  to  be  adapted  in  an  all-zero  form. 
The  reformulation  restores  the  quadratic  nature  of  the  minimization  problem  to 
that  encountered  in  the  all-zero  case  and  allows  easy  adaptation  using  LMS  [3], 
RLS  (24],  sample  matrix  inversion  [10]  or  any  other  least-squares  algorithm.  A 
least-squares  adaptive  lattice  implementation  of  the  equation-error  array  method 
has  been  worked  out  in  [25], 

Refer  to  Fig.  11  and  assume  for  the  moment  that  the  linear  filters  are  time 
invariant.  The  transfer  function  from  point  d  to  the  output  is  unity  since  the 
all-pole  and  all-zero  filters  cancel.*  This  cancellation,  combined  with  the  fact 
that  the  look-direction  signal  does  not  appear  in  the  zk  inputs,  allows  the  look- 
direction  signal  to  pass  through  the  array  undistorted.  The  transfer  function 

Bk(z) 

from  point  zk  to  the  output  is  - -  and  contains  both  poles  and  zeros. 

1  +  z~'A(z) 

Instead  of  minimizing  the  overall  array  output  (which  contains  the  all-pole  filter), 
our  approach  minimizes  the  signal  appearing  at  the  input  to  the  all-pole  filter.  In 
the  system  identification  literature,  this  signal  is  often  referred  to  as  the  equation 
error  [23].  Justification  for  minimizing  the  equation  error  rather  than  the  true 
output  error  is  that  in  many  cases  by  making  the  equation  error  small,  the  out¬ 
put  error  is  also  made  small.  Further  discussion  of  the  differences  between  equa¬ 
tion  error  and  output  error  can  be  found  in  [27], 

The  most  important  aspect  of  the  equation-error  approach  is  that  the  mean 

square  of  the  equation  error  is  a  quadratic  function  of  the  pole  polynomial 

*Due  to  their  adaptive  nature,  the  array  filters  are  time  varying.  It  is  assumed,  howev¬ 
er,  that  the  adaptation  time  constant  is  long  compared  to  the  time  constant  of  the  filters. 

In  this  case,  the  array  filters  can  be  considered  time  invariant. 


Figure  11.  A  pole-zero  adaptive  array  processor  with  unity  gain  look  direction  constraint 


k 
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coefficients.  The  convergence  rate  of  the  system  can  be  made  as  fast  as  the 
signal-to-noise  ratio  will  allow  using  sample  matrix  inversion  of  other  exact- 
least-squares  algorithms.  This  is  not  the  case  in  conventional  pole-zero  adaptive 
algorithms,  such  as  RNIL,  which  minimize  the  output  error.  Such  algorithms 
usually  exhibit  slow  convergence  and  cannot  guarantee  globality  of  the  con¬ 
verged  solution. 

To  convert  the  pole-zero  array  of  Fig.  11  into  an  equivalent  array  which  uses 
unconstrained  adaptation,  observe  that  the  transfer  function  from  sensor  k  to  the 
output  is  given  by 


Wk(z) 


^  [i+z-M0)|-Bi(0+  fl1+1(0 

1+  2_1A(z) 


(52) 


where  fl0(c)  =  )  =  0.  It  is  easily  seen  that  eq.  (52)  can  be  equivalently 

written  as 


Wk{*) 


nk0  +  z  XNk  (*) 
1+  z~^  A.  (z) 


(53) 


m 

where  n*o  =  1  and  the  Nk  (z)  polynomials  are  unconstrained.  This  leads  to 

the  alternate  structural  implementation  shown  in  Fig.  12.  The  similarity 
between  this  system  and  the  Frost  system  shown  in  Fig.  1  is  apparent.  Instead 
of  constraining  the  sum  of  all  the  filters  to  be  unity,  only  the  sum  of  the  first 
coefficients  has  been  constrained.  A  look-direction  signal  will  experience  a 
transfer  function  from  input  to  the  adaptation  error  point  given  by 


Figure  12.  Alternate  realization  of  a  pole-zero  array  processor 
with  unity-gain  look-direction  constraint. 
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I  +  »'*  S  •%’<«)  •  (54) 

This  transfer  function  is  inverted  by  the  all-pole  output  filter  to  generate  the 
array  output  and  restore  the  unity- gain  look-direction  constraint. 

It  should  also  be  mentioned  that  the  array  processing  structure  shown  in 
Fig.  12  could  be  used  in  non-adaptive  applications.  For  instance,  it  could  be  use¬ 
ful  in  an  acoustic  application  where  the  objective  was  to  listen  in  one  direction 
while  placing  deep  wideband  nulls  in  other  directions.  Knowledge  of  the  direc¬ 
tions  and  spectra  of  the  sources  could  be  used  to  calculate  the  coefficients  of  the 
array.  These  coefficients  would  then  be  built  into  a  fixed  non-adaptive  pole-zero 
array  processor. 

Stability  of  the  inverse  filter 

One  issue  yet  to  be  addressed  is  that  of  stability  of  the  inverse  filter.  In 
some  situations,  it  is  possible  that  the  roots  of  the  polynomial  associated  with  the 
inverse  filter  can  move  outside  the  unit  circle.  When  this  occurs,  the  inverse 
filter  as  shown  in  Fig.  11  will  become  unstable  and  its  output  will  begin  to  grow 
without  bound. 

A  pole  polynomial  with  roots  both  inside  and  outside  the  unit  circle 
represents  a  stable  filter  only  in  a  noncausal  sense.  It  has  an  impulse  response 
which  is  two-sided,  extending  from  k=- oo  to  k= oo.  Such  a  response  can  be 
approximately  realized  in  a  delayed  form.  For  instance,  assume  that  the  anti- 
causal  part  of  the  above  impulse  response  extends  only  to  k=-N  and  is  zero 
from  k=-N  to  k=-oo.  By  introducing  an  N-stcp  delay,  this  impulse  response 
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can  be  time-shifted  so  that  it  is  entirely  causal  and  thus  realizable.  In  general, 
the  transfer  function  can  be  factored  into  the  product  of  a  purely  causal  part 
with  all  poles  inside  the  unit  circle  and  a  purely  anticausal  part  with  all  poles 
outside  the  unit  circle.  The  causal  part  can  be  exactly  realized.  The  anticausal 
part,  however,  has  an  impulse  response  that  theoretically  extends  back  to  the  minus 
infinity.  Practically,  the  anticausal  inpulse  response  has  a  finite  time  duration 
whose  length  is  determined  by  the  proximity  of  the  transfer  function  poles  to  the 
unit  circle.  The  closer  the  poles  are  to  the  unit  circle,  the  longer  the  impulse 
response. 

One  method  for  overcoming  the  possible  instability  of  the  inverse  filter  is  to 
approximate  it  using  an  all-zero  filter.  In  order  to  yield  an  accurate  approxima¬ 
tion,  the  all-zero  filter  needs  to  be  of  high  order  and  must  be  approximated  in  a 
delayed  form.  The  consequences  of  introducing  the  delay  are  minimal.  The 
consequences  of  using  a  high  order  all-zero  filter  as  an  approximate  inverse  are  a 
bit  more  severe.  Additional  hardware  is  required  to  both  form  and  implement 
the  all-zero  approximation.  At  first  glance,  it  might  seem  as  though  making  a 
high-order  all-zero  approximation  defeats  the  original  purpose  of  using  poles  to 
economize  on  weights.  However,  note  that  the  adaptive  algorithm  still  operates 
on  a  reduced  number  of  coefficients.  Formation  of  the  approximate  inverse  can 
be  done  off-line  and  does  not  introduce  complexity  into  the  adaptation.  The 
speed  of  convergence  of  the  overall  system  will  not  be  affected.  Details  of  form¬ 
ing  an  all-zero  inverse  approximation  are  given  in  reference  [26]. 

The  preferred  method  of  overcoming  the  inverse  filter  instability  is  to 


prevent  the  pole  polynomial  roots  from  ever  going  outside  the  unit  circle.  This 
can  be  done  in  one  of  two  ways.  The  first  way  is  to  constrain  the  polynomial  to 
remain  minimum  phase  (i.e.,  keep  all  roots  within  the  unit  circle)  during  adapta¬ 
tion  using  a  constrained  adaptive  algorithm  which  monitors  the  polynomial  roots. 
The  second  way  is  to  introduce  a  bulk  delay  into  the  desired  response  path  of 
Fig.  11.  It  turns  out  that  increasing  the  bulk  delay  causes  the  roots  of  the  con¬ 
verged  pole  polynomial  to  move  closer  to  the  origin.  Choosing  a  delay  too  large 
causes  the  resulting  solution  to  lose  quality  and  is  not  desirable.  Choosing  a 
delay  too  small  results  in  a  noncausai  solution  and  an  unstable  inverse  filter.  The 
idea  is  thus  to  choose  the  smallest  delay  which  causes  all  the  roots  to  lie  within 
the  unit  circle. 

Currently,  we  are  experimenting  with  several  on-line  methods  for  adaptively 
choosing  the  above  bulk  delay.  A  brute  force  method  is  to  have  a  series  of 
almost  identical  systems  all  operating  on  the  same  input  data  but  each  using  a 
different  desired  response  delay.  A  simple  criterion  is  then  used  to  choose  the 
stable  system  having  the  smallest  delay.  Another  method  is  to  use  a  single  sys¬ 
tem  having  a  variable  delay  that  can  be  incremented  or  decremented  depending 
on  the  stability  of  the  pole  polynomial. 

Look-direction  signal  bias 

The  steady-state  Wiener  solution  of  the  conventional  Frost  array  is  not 
affected  by  the  presence  of  a  look-direction  signal.  This  is  readily  apparent  by 
referring  to  Fig.  6  and  noting  that  the  x-input  does  not  contain  look-direction 
signal.  Therefore  the  look-direction  signal  cannot  contribute  to  either  the  cross- 
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correlation  between  the  desired-response  input  d  and  the  x-input  or  the  auto¬ 
correlation  of  the  x-input. 

The  Wiener  solution  of  the  pole-zero  Frost  array  is,  however,  affected  by  a 
look-direction  signal.  To  demonstrate  this  effect,  refer  to  Fig.  11  and  assume 
that  the  only  source  of  signal  originates  from  the  look  direction.  In  this  condi¬ 
tion,  the  x-input  will  be  zero  and  the  d-input  will  consist  of  look-direction  alone. 
During  adaptation,  the  1+  z^A[z)  filter  will  act  as  a  linear  predictor  and  try  to 
reduce  the  look-direction  signal  power.  If  the  signal  were  narrowband,  the  pred¬ 
ictor  would  in  fact  form  a  notch  filter  at  the  signal  frequency.  Distortion  of  the 
look-direction  signal  would  not  occur  since  the  notch  filter  response  would  be 
inverted  by  the  inverse  output  filter  to  restore  the  look-direction  signal  at  the 
array  output.  The  main  concern  is  not  signal  distortion,  but  the  loss  of  degrees 
of  freedom  caused  by  attempting  to  internally  null  the  look-direction  signal.  This 
degradation  is  referred  to  as  look-direction  signal  bias. 

In  situations  where  adaptation  can  be  controlled  so  that  it  occurs  only  when 
the  look-direction  signal  is  not  present,  look-direction  signal  bias  will  not  occur. 
This  is  the  case  in  many  radar  or  spread-spectrum  applications  where  the  desired 
signal  is  pulsed  or  sequenced  in  a  manner  which  is  known  to  the  receiver.  Also, 
look-direction  signal  bias  will  not  cause  a  significant  degradation  in  the  array  per¬ 
formance  if  the  look-direction  signal  power  is  substantially  less  th  m  the  interfer¬ 
ence  power.  In  a  system  where  the  desired  signal  is  continuous  and  of  significant 
power,  a  solution  to  the  bias  problem  must  be  found. 

One  method  of  alleviating  look-direction  signal  bias  is  shown  in  Fig.  13.  The 
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method  is  based  on  an  idea  proposed  by  Duvall  for  eliminating  dynamic  signal 
cancellation  [20,21].  For  simplicity,  only  a  three-element  example  has  been 
shown.  The  basic  idea  is  to  process  the  set  of  signals  from  the  sensors  in  a 
manner  which  eliminates  the  look-direction  signal  but  which  keeps  the  phase  rela¬ 
tionships  among  the  new  set  of  signals  the  same  as  that  prior  to  processing.  The 
processed  set  of  signals  is  fed  into  an  adaptive  “master”  array  where  the  Wiener 
solution  cannot  be  affected  by  the  look-direction  signal.  The  weights  of  the  mas¬ 
ter  array  are  copied  into  the  “slave”  array  which  operates  on  the  original  signal 
set.  Because  of  the  preservation  of  phase,  nulls  formed  using  the  slaved  array 
will  be  in  the  same  directions  as  those  formed  in  the  master.  For  more  details  see 
reference  [20], 

To  analyze  the  effect  of  preprocessing  on  the  Wiener  solution,  we  define  an 
effective  interference  environment  that  incorporates  the  preprocessor  into  the 
interference  sources.  Each  interference  source  experiences  an  angle-dependent 
frequency  transformation  when  passed  through  the  preprocessor.  The  frequency 
response  of  this  transformation  is  determined  by  the  subtraction  of  signals 
appearing  at  adjacent  sensors.  Again  we  assume  that  the  array  of  sensors  is 
linear  and  uniformly  spaced.  A  directional  interference  source  will  experience  a 
delay,  A,  between  two  adjacent  sensors  that  is  determined  by  the  interference 
bearing  angle.  The  frequency  response  of  the  top  subtractor  in  the  preprocessor 
is  given  by 

//(wA)  =  1  -  t>u£k  .  (55) 


Proceeding  down  the  array,  each  subtractor  exhibits  this  same  response  plus  an 


additional  delay.  The  idea  is  to  use  the  //( lj)  transformation  to  modify  the 
interference  source  and  conceptually  eliminate  the  preprocessor  from  the  array. 


Notice  that  as  the  interference  angle  approaches  the  look  direction,  the  mag¬ 
nitude  of  the  preprocessor  transfer  function  goes  to  zero.  The  look-direction  sig¬ 
nal  will  be  removed  and  an  interference  source  appearing  off  the  look  direction 
will  be  somewhat  attenuated  and  high-pass  filtered.  In  this  manner,  the  prepro¬ 
cessor  affects  the  Wiener  solution  by  slightly  modifying  the  weighting  applied  to 
the  interference  source. 

Sim  illations 

Figure  14  compares  the  frequency  response  of  the  optimal  weight  at  each 
sensor  with  the  frequency  response  of  the  converged  array  filters  obtained  using 
the  proposed  pole-zero  adaptive  array  processor.  In  this  simulation,  the  look- 
direction  signal  was  assumed  to  be  time-sequenced  in  a  manner  which  allowed  the 
array  to  adapt  only  when  the  signal  was  not  present.  Thus,  look-direction  signal 
bias  was  not  a  concern.  As  in  the  previous  simulation,  the  interference  source 
had  a  normalized  center  frequency  of  7r/2,  a  100%  relative  bandwidth,  and  a 
bearing  angle  45  degrees  off  the  look  direction.  Again  there  was  no  thermal  noise. 
The  array  filters  had  2  zeros  each  and  2  common  poles.  This  translates  into  a 
total  of  8  filter  coefficients  compared  to  the  10  coefficients  used  in  the  Frost  array 
simulation.  Notice  that  the  pole-zero  filter  response  plotted  in  Fig.  14  is  a  much 
closer  fit  to  the  optimal  response  than  the  conventional  all-zero  filter  response 
plotted  in  Fig.  7. 


Comparison  of  frequency  response  of  optimal  array  weights  (cross-hatched  marks) 
with  a  2  pole  -  2  zero  approximation,  (a)  Sensor  1  (b)  Sensor  2. 
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For  further  comparison,  Fig.  15  shows  the  broadband  antenna  pattern  and 
the  frequency  response  of  the  array  in  the  signal  and  interference  directions.  The 
interference  null  is  about  40  dB  below  the  look-direction  gain.  This  represents  a 
20  dB  improvement  over  a  conventional  Frost  system  having  essentially  the  same 
number  of  variable  coefficients. 

Fig.  16  shows  the  nulling  capabilities  of  a  pole-zero  Frost  array  processor  as 
a  function  of  the  interference  bearing  angle.  The  conditions  of  this  experiment 
were  the  same  as  those  given  above  except  that  -30  dB  of  white  thermal  noise 
was  added  to  each  sensor.  Comparing  this  graph  with  the  one  presented  in  Fig. 
9,  we  see  that  throughout  the  field  of  view,  the  pole-zero  Frost  array  outperforms 
the  conventional  Frost  array.  Notice  that  as  the  interference  bearing  angle 
approaches  the  look  direction,  the  null  depth  rises  to  0  dB  and  the  unity  gain 
look-direction  constraint  is  preserved. 

A  final  set  of  simulations,  shown  in  Fig.  17,  compares  the  converged  broad¬ 
band  antenna  pattern  of  a  conventional  Frost  array  with  the  pattern  of  a  pole- 
zero  Frost  array.  Both  arrays  consisted  of  three  colinear  uniformly  spaced  sen¬ 
sors.  Incident  on  both  arrays  were  two  mutually  uncorrelated  wideband  interfer¬ 
ence  sources  each  having  a  normalized  center  frequency  of  ir/2  and  a  100%  rela¬ 
tive  bandwidth.  The  conventional  Frost  array  had  6  zeros  per  sensor.  The  pole- 
zero  Frost  array  had  4  zeros  per  sensor  and  4  poles  common  to  all  sensors.  Thus, 
both  systems  had  essentially  the  same  number  of  variable  coefficients.  Using  an 
exact  least  square  algorithm,  an  equal  number  of  coefficients  translates  into 
approximately  equal  convergence  rates.  Thus  the  pole-zero  Frost  antenna  pat- 
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Broadband  thermal  noise  power  and  directional  interference  power  at  output  of 
a  two-element  pole-zero  Frost  array  plotted  as  a  function  of  intersensor  delay. 


-  55- 


tern  shown  in  Fig.  17b  required  the  same  amount  of  data  to  reaoh  convergence  as 
did  the  conventional  Frost  pattern  shown  in  Fig.  17a. 
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VI.  FUTURE  RESEARCH 

The  work  presented  in  this  paper  has  demonstrated  the  remarkable  utility  of 
the  proposed  pole-zero  adaptive  array  processor.  As  with  any  important  new 
research  result,  this  paper  has  opened  up  a  wide  area  for  future  research.  Listed 
below  are  a  few  topics  in  need  of  further  study. 

The  application  of  the  pole-zero  adaptive  processor  to  nonconformal  arrays 
should  be  analyzed  and  simulated.  More  work  should  also  be  conducted  concern¬ 
ing  stabilization  of  a  nonminimum-phase  pole  polynomial.  The  covariance  matrix 
eigenvalue  spread  should  be  analyzed  in  order  to  determine  convergence  proper¬ 
ties  of  the  pole-zero  adaptive  array  when  implemented  using  the  LMS  algorithm. 
Also,  the  possibility  of  improved  convergence  when  implementing  the  pole-zero 
processor  with  an  adaptive  lattice  structure  should  be  explored. 
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vn.  CONCLUSIONS 

In  this  paper,  a  new  wideband  adaptive  array  processing  structure  was 
presented.  Using  a  derivation  of  the  optimal  wideband  array  weighting  function, 
it  was  conjectured  that  filters  possessing  both  poles  and  zeros  would  be  capable  of 
improving  the  performance  of  the  nulling  array.  The  “equation-error”  approach 
was  developed  and  used  as  a  method  for  adapting  the  poles  and  zeros  of  the 
array  filters.  The  causes  and  cures  of  two  undesirable  attributes  of  the  equation- 
error  approach  were  discussed.  Simulations  like  the  one  presented  in  Fig.  17 
demonstrate  that  the  new  pole-zero  array  processor  can  dramatically  improve  the 
wideband  nulling  capability  of  the  array,  providing  sharper  and  deeper  nulls 
while  using  the  same  total  number  of  weights. 
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